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Abstract 

In this letter, we extend the known results for the QCD potential between a static 
quark and its antiquark by computing the two-loop corrections to the colour-octet 



X 

. p. , state. 

PACS numbers: 12.38.Bx, 12.38.-t 

The QCD potential between a static quark and its antiquark has for a long time been 
used as a probe of the fundamental properties of the strong interactions such as asymp- 
totic freedom and confinement p. Historically, the potential for a quark-antiquark pair 
in the colour-singlet state attracted the most attention because it is a basic ingredient 
in the theory of heavy quarkonium and, therefore, of primary phenomenological interest. 
Nowadays, however, there is growing interest in its colour-octet counterpart. The latter 
naturally appears in effective-theory calculations of high-order corrections to the heavy- 
quarkonium spectrum and decay rates through the so-called ultrasoft contribution j2]. 
Moreover, it determines the properties of glueballinos and is necessary for the analysis of 
gluino-antigluino threshold production |3|4j . It is also used in lattice QCD for studying 
the behavior of strong interactions at long distances and the interplay between pertur- 
bative and non-perturbative physics [3]. This requires knowledge of the corresponding 
perturbative corrections which, in contrast to the colour-singlet not available 



beyond one loop. In the present letter, we fill this gap and compute the C(a^) correction 
to the colour-octet static potential. 

The perturbative expansion of the colour-singlet potential reads 
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where the first term corresponds to the Coulomb potential. The one-loop coefficient, 
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has been known for a long time |5|6j . while the two-loop coefficient, 02, has only recently 
been found (7i|i8ip9j . In Ref. 0, the result of Ref. jBj was confirmed. 
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where C, is Riemann's zeta function, with value C(3) = 1.202057.... Here, Ca = N 
and Cp = (iV^ — 1)/(2A^) are the eigenvalues of the quadratic Casimir operators of the 
adjoint and fundamental representations of the SU{N) colour gauge group, respectively, 
Tp = 1/2 is the index of the fundamental representation, and ni is the number of light- 
quark flavours. The modified minimal-subtraction (MS) scheme for the renormalization 
of as is implied. The logarithmic term of 0{al) in Eq. (0) refiects the infrared divergence 
of the static potential The particular form of the logarithmic term corresponds to 
dimensional regularization ^T]. The corresponding infrared-divergent term is cancelled 
against the ultraviolet-divergent one of the ultra-soft contribution [2^ in the calculation of 
the physical heavy-quarkonium spectrum |ll|12j . The non- logarithmic third-order term, 
03, is still unknown. 

The perturbative expansion of the potential for the colour-octet state can be cast in 
the form 
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where the one-loop coefficient is the same as in the colour-singlet 
two- loop coefficient, however, differs by a finite renormalization-independent term, 
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Our result is 
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= (tt^ - 127r2) C\ + 0(d - 4) , (6) 

where d is the space-time dimension. The non-logarithmic part of the three-loop coeffi- 
cient, Og, is still unknown. It is instructive to look at the numerical size of the corrections. 
For = 3 one obtains ba, = —189.2. At the same time, we have a, = 155.8(211.1, 268.8) 
and ai = 4.778(5.889, 7.000) for rii = 5(4, 3). Thus, in the colour-octet case, the two-loop 
correction is significantly smaller than for the colour-singlet configuration. Depending on 
rii, it even changes sign. 

In the remaining part of this letter, we wish to describe two independent ways that 
have been used to evaluate 6a2- The ffist method proceeds along the lines of the analysis 
j9|llj based on the threshold expansion jT^. In general, the threshold expansion is the 
proper framework for performing calculations involving a heavy quark- ant iquark system. 
It provides rigorous power-counting rules and natural definitions of the formal expressions 
obtained in the perturbative analysis of the nonrelativistic effective theory. The correc- 
tions to the static potential only arise from the soft regions of the loop integrals, which 
are characterized by the following scaling of the loop momenta: Iq ~ \l\ ~ \q\. Thus, the 
calculation of the coefficients and a° can be performed in the static limit of NRQCD, 
niq — > oo. 

Due to the exponentiation of the static potential jSj, the coefficients of the colour- 
singlet state only receive contributions from the maximally non-Abelian parts, leaving 
aside the terms involving n/. The selection of these parts effectively retains the contri- 
butions of the soft region, as the appearance of the Abelian colour factor Cp indicates 
the presence of a Coulomb pinch and thus implies that at least one loop momentum is 
potential. The latter contributions just represent iterations of the lower-order potential 
and, therefore, should be excluded from the potential itself. In the nonrelativistic effec- 
tive theory, these iterations are taken into account in the perturbative solution of the 
Schrodinger equation about the Coulomb approximation. These contributions refer to 
dynamical rather than static heavy-quark and -antiquark fields, and the Coulomb pinch 
singularities we encounter in the static-limit calculations are resolved by keeping a finite 
mass in the nonrelativistic heavy-quark propagator. 

The analysis of the colour-octet state is more involved, since, in this case, the Coulomb 
pinches come with all possible colour factors and cannot be removed by selecting the 
maximum non-Abelian ones. Thus, the separation of the Coulomb pinches should be 
performed explicitly. They appear in the Feynman diagrams involving the product of the 
nonrelativistic quark and antiquark propagators, 
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In this case, after expanding the quark propagator in l/m^, one obtains ill-defined prod- 
ucts like 



Thus, separating the soft and potential regions is unavoidable.^ In the soft region, the 
pole contributions of the quark and antiquark propagators have to be excluded, and the 
product in Eq. (jSI) should actually be defined to be its principal value. 
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In the potential region, the quark and antiquark propagator poles produce contributions 
of the form 
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where the 1/v Coulomb singularity shows up explicitly. After integration over ko, Eq. (|T0|) 
yields the nonrelativistic Green function of the free Schrodinger equation. Only Eq. 
should be taken into account in the calculation of the static potential. 

At one-loop, there is only one diagram involving a Coulomb pinch, namely, the planar 
box, which has the colour factor Cp for the colour-singlet state. Picking up the soft 
contribution, i.e., using the principal- value prescription of Eq. Q to define Eq. (jHl), we find 
the planar box to cancel the Cp part of the non-planar box, which in total is proportional 
to Cp — CfCa/2. This explicitly demonstrates the exponentiation of the one-loop colour- 
singlet static potential in momentum space. However, we can also turn things around and 
express the planar box with Coulomb pinches through the well-defined non-planar box by 
actually requiring the cancellation of the Cp terms in the sum of all one-loop diagrams, 
as is dictated by the exponentiation. The result for a° as given above is then obtained by 
simply replacing the colour- singlet colour factor by the colour-octet one. 

This strategy carries over to two loops. Here, we have diagrams with zero, one, or 
two Coulomb pinches. For the diagrams without Coulomb pinch, the contribution to 
02 is obtained by adopting the correct colour factor. We divide the Feynman diagrams 
with Coulomb pinches into those that have two quark and two antiquark propagators 
(cf. Fig. and the rest. The latter ones are treated directly using the principal-value 
prescription of Eq. 0. For the former, however, it is simpler to use the exponentiation, 
which requires that the diagrams contributing to the colour factors CaCp and Cp sum up 
to zero in the colour- singlet case. This leads to two equations for the diagrams suffering 
from Coulomb pinches, namely, those shown in Figs. ^ and b, which can be solved. This 
provides a result in terms of the diagrams in Figs.QJ^ and d, which are free of pinches. After 

^Note that, for the diagrams without Coulomb pinches, the separation of the soft and potential regions 
is ambiguous and even gauge dependent. In such diagrams, the nonrelativistic quark and antiquark 
propagators can be safely expanded in 1/mq. 
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Figure 1: Two- loop Feynman diagrams with ((a) and (b)) and without ((c), (d), and (e)) 
Coulomb pinches that contribute to 5a2- 



adopting the colour factors corresponding to the colour-octet configuration, one obtains 
the contributions to the results given in Eqs. (0) and We wish to mention that the 
calculation was performed in the general covariant gauge and that the dependence on the 
gauge parameter was found to cancel out in the final result. 

The second method to compute V° proceeds along the lines of Ref . Pj . While in the 
above, we had to assume exponentiation of the colour- singlet potential, we will now relax 
that assumption. The reason is that, although exponentiation is plausible to all orders of 
perturbation theory, the proof given in Ref. [5] holds for the singlet potential in Abelian 
theories only. 

As a starting point, we now expand the logarithm of the (T x R) Wilson loop spanned 
by the static quark- ant iquark pair at distance R through 0{al). Taking the limit T —>■ oo 
(which, in a diagrammatic sense, 'cuts' the Wilson loop twice and restores translational in- 
variance in the temporal direction, hence guaranteeing energy conservation at the vertices 
and leading to simple momentum-space Feynman rules) and inserting SU {N) generators 

into the purely spatial Wilson lines to obtain the colour-octet potential to this or- 
der (for a manifestly gauge-invariant definition, see Ref. ^1]), we now explicitly keep 
disconnected as well as one-particle-reducible diagrams in our expansion. 

At this point, the general structure of the expansion involves (products of) up to 
two-loop four-point functions of static quarks (cf. Fig. P). After Fourier transforming to 
momentum space, we can choose a special point to evaluate these four-point functions, 
since the potential, of course, only knows about the distance R of the qq pair, which in a 



5 



momentum space representation translates into the momentum transfer \q\ between the 
upper and lower lines in Fig. ^ Hence, effectively, we have to compute two-point functions 
with external static quarks, external momentum q = {0,q), and internal static quarks, 
gluons, ghosts and light quarks, with the additional occurrence of a static (anti-) quark - 
gluon two-point vertex, resulting from the special kinematics. 

After performing the colour algebra and exploiting symmetries of the integrals occur- 
ring in the expansion, all integrals which might give rise to pinch singularities, and had 
to be treated with caution in our first approach, cancel exactly. Thus, we are left with 
the task of computing a class of two-loop two-point integrals for which there exists a 
generic algorithm ^8ifl5] . based on integration by parts (IBP) ^HI- The implementation 
in Ref . [T3] (see also Chapter 6 of Ref . [T7] ) is based on Ref . [TH] . 

Having generated the relevant set of diagrams and reduced the occurring Feynman 
integrals to the set of two-point functions described above, we now employ the reduction 
algorithm, which maps them to a (small) set of so-called master integrals, multiplied by 
rational functions in the dimension d. At this stage, we observe cancellation of the gauge- 
parameter dependence, serving as a check for the reduction. As an additional strong 
check, we use our implementation [1^] of the strategy to solve a truncated set of IBP 
relations, based on lexicographic ordering of integrals [2U]. 

The set of (massless, two-point) master integrals is known analytically in terms of 
Gamma functions, for generic dimension d, as given, e.g., in Ref. |H]. Expanding prefactors 
as well as master integrals about d = 4 — 2e and renormalizing the gauge coupling, we 
again arrive at Eq. ©. 

To conclude, we have evaluated the C(a^) correction to the colour-octet static po- 
tential using two independent techniques. Both evaluations are in agreement, giving us 
confidence in our main result, Eq. 
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